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Аннотация 

Введение. Производительность и надежность работы героторных гидромашин зависит от геометрических 
параметров профиля циклоидального зацепления. Существующие методики расчета и оптимизации параметров 
профиля громоздкие, многокритериальные, сложные для практического применения. Поэтому актуальной 
представляется проблема создания методики расчета параметров профиля рабочего органа героторной машины, 
пригодной для инженерных расчетов на этапе эскизного проектирования. В связи с этим целью данной работы 
являлась модернизация методики проектирования геометрии профилей торцевого сечения гипоциклоидальных 
зубчатых зацеплений, используемых в героторных гидравлических машинах, и анализ возможностей их 
оптимизации при предварительном проектировании. В ходе исследования была использована система 
компьютерной математики Ма#са4, проведены численные эксперименты для изучения влияния 
геометрических параметров профилей на производительность и работоспособность героторной гидравлической 
машины. На основе полученных и проанализированных данных выработаны рекомендации по проектированию 
оптимальных профилей торцевого сечения героторных гидравлических машин. 

Материалы и методы. Материалы включают в себя известные методики расчета параметров профиля, 
основанные на применении классических формул эквидистант гипоциклоид, используемых для очерчивания 
профилей зубьев рабочих органов героторных машин. Основной метод исследования — моделирование 
профиля героторной машины с помощью системы компьютерной математики Мафса4. Получены расчетные 
данные для выбранных диапазонов варьируемых параметров, обработанные методом регрессионного 
однофакторного анализа. 

Результаты исследования. Разработан алгоритм анализа гладкости профилей зубьев. Определены два целевых 
параметра: площадь сечения торцевого профиля, влияющая на производительность, и наименьший 
приведенный радиус контакта, определяющий работоспособность рабочего органа. Предложена методика по 
расчету целевых параметров на ранней стадии проектирования. Получен ряд оптимальных значений 
параметров профиля по критериям производительности и работоспособности героторной машины. Построены 
зависимости, позволяющие установить оптимальные значения параметров профиля на стадии проектирования. 
Обсуждение и заключение. Разработанная авторами методика дает возможность на этапе проектирования 
рабочего органа получить оценку производительности и работоспособности героторной гидромашины. 
Результаты исследований могут быть использованы в машиностроении при проектировании героторных 
гидромашин для улучшения их технико-эксплуатационных характеристик. 
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Киреев С.0О. и др. Оптимизация геометрических характеристик циклоидальных профилей героторных гидромашин 


Введение. Героторные гидравлические машины все чаще стали применяться в строительстве и при 
эксплуатации нефтяных и газовых скважин, причина — их конструктивно-эксплуатационные преимущества, 
основными из которых являются низкая неравномерность подачи для насосов или крутящего момента для 
двигателей, компактные габаритные размеры. 

В научной литературе большое место занимают исследования, посвященные работе героторных 
гидромашин. Сильное влияние геометрических характеристик профиля ротора на производительность и 
надежность винтовых героторных гидромашин отмечают Д. и Ф. Балденко [1]. Основным узлом такой машины 
является многозаходный винтовой героторный механизм — цилиндрическая планетарная зубчатая передача 
внутреннего зацепления, состоящая из статора и ротора, с разницей в числах зубьев, равной единице. Профиль 
статора в торцовом сечении контактирует с профилем ротора, образуя замкнутые площадки. Авторы 
проанализировали также технологии изготовления зубчатых передач. 

Р. Айер, В. Джатти исследовали эффективность применения математических моделей при конструировании 
винтовых насосов [2]. Они представили также результаты моделирования различных типоразмеров 
двухвинтовых насосов при перекачивании воды, смеси воды и воздуха. Суан-Фыонг Данг и Хонг-Сок Парк 
предложили процесс оптимизации САП-модели с использованием результатов моделирования САЕ и вычислений 
в Ма@саа [3]. Другие авторы (А. Лебедев, С. Киреев, М. Корчагина, С. Власкин, А. Ефимов) в своих публикациях 
отметили эффективность проектирования разнообразных объектов и процессов на основе параметризации и 
последующей оптимизации конструкции по одному или нескольким целевым параметрам [4—6]. 

В совместной работе Ф.Д. Балденко и Яо Ян подчеркнули, что площадь живого сечения не зависит от 
относительного положения профилей ротора и статора и определяет производительность насоса, 
проанализировали влияние безразмерных геометрических параметров профиля на площадь живого сечения [7]. 

Лэй Чжэн, Сяодун Ву, Гоцин Хань, Хуачан Ли, Йи Цзо и Дэйк Чжоу разработали конечно-элементную 
модель рабочего органа винтового насоса с поверхностным приводом для нефтегазовой промышленности [8]. 
По результатам моделирования и экспериментальных проверок ими был сделан вывод, что зазоры в 
сопряжении и толщина статора являются двумя основными факторами, влияющими на оптимизируемый 
целевой параметр — объемную эффективность насоса. 

Новые данные получены Х. Гамбоа, А. Оливет, С. Эспин при моделировании рабочего органа 
однозаходного винтового насоса [9]. Эти авторы определили существенную зависимость между площадью 
зазора проскальзывания и перепадом давления, связанным с механическими свойствами материала статора. 

Яо Ян, Ф. Балденко исследовали влияние безразмерных коэффициентов торцового профиля героторной 
гидромашины на характеристики циклоидного зацепления, что позволяет проводить оптимизацию на основе 
геометрических и кинематических критериев [10]. И.А. Лягов, Ф.Д. Балденко, А.В. Лягов, В.У. Ямалиев, 
А.А. Лягова предложили методику выбора оптимальной конфигурации силовых секций при проектировании 
секционных забойных двигателей, отметив актуальность проблемы уменьшения габаритных параметров длины 
винтового насоса [11]. 

Ю.А. Коротаев, А.Н. Алпатов, А.В. Соболев и др. пришли к выводу, что наиболее эффективным является 
профилирование зацепления от исходного контура рейки, очерченного эквидистантой укороченной 
циклоиды [12]. В отношении профилей, полученных от эквидистанты обыкновенной гипоциклоиды, авторы 
указали на их основной недостаток: они строятся как специальные и не предусматривают натяг в зацеплении. 
Ю.А. Коротаев и Д.А. Голдобин обратили внимание на проблему, связанную с необходимостью повышения 
производительности винтовых насосов, и возникающие при этом недопустимые условия жесткости при 
производстве статоров и роторов [13]. 

А.Ф. Миникаев, В.А. Пронин, Д.В. Жигновская, Ю.Л. Кузнецов, Ф.Д. Балденко, А.Е. Ковалёнок разработали 
методику исследования напряженно-деформированного состояния статора при контактном взаимодействии с 
ротором винтового рабочего органа гидромашины с применением САЕ-моделирования и утверждают, что 
классическая теория Герца о контактном взаимодействии не может быть применена для определения 
контактных напряжений в роторе и статоре из-за действия ряда особенностей: изменения кривизны 
соприкасающихся поверхностей ротора и статора, различной толщины эластичной обкладки по длине 
окружности корпуса, наличия предварительного натяга в паре, сложного планетарного движения ротора, 
сопровождающегося сочетанием трения качения и скольжения, смещения вектора радиальной силы 
относительно нормали в точке контакта. 

Подводя итог обзору научных публикаций на тему проектирования героторной гидромашины, необходимо 
отметить, что большинство исследователей в данной области стремятся повысить эффективность работы этой 
машины, увеличивая площадь живого сечения и сокращая габаритные размеры рабочего органа. Достижение 
оптимальных результатов происходит за счет применения и варьирования коэффициентов внецентроидности, 
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формы зуба, относительной скорости скольжения, смещения рейки ит. д. Методики расчета громоздкие, их 
сложно применять на практике. Торцевые профили строятся, как правило, на основе укороченных гипоциклоид, 
что приводит к возникновению проблемы с обеспечением кинематической точности движения, контактного 
взаимодействия ротора и статора. Некоторые авторы (Ю.А. Коротаев, А.Н. Алпатов, А.В. Соболев и др.) указали 
на невозможность создания натяга в случае применения обыкновенной гипоциклоиды, но натяг может быть 
создан за счет изменения эквидистантного расстояния профилей ротора и статора. 

Материалы и методы. Для получения профилей ротора и статора в торцевом сечении рабочего органа 
героторной гидравлической машины зададимся следующим набором исходных параметров (рис. 1): 

1) Е, — контурный радиус профиля, определяющий диаметр окружности С!, описывающей профиль 
статора (мм). Фактически этот параметр задает радиальный габарит героторной машины, что особенно важно в 
нефтегазовой промышленности, когда проектируемый объект необходимо разместить, например, в стволе 
эксплуатационной колонны скважины; 

2) =, — число зубьев ротора; 

3) В. — радиус ролика, определяющий расстояние эквидистантной линии профилей ротора и статора 
относительно соответствующих гипоциклоид (мм); 

4) А — натяг в сопряжении ротора и статора (мм). 

Напомним принцип образования профилей гипоциклоидального зацепления, который заключается в том, 
что внутри окружности С, (основная окружность гипоциклоиды статора) перекатывается без скольжения 


окружность С, (основная окружность гипоциклоиды ротора). Центр окружности С, всегда находится на 


расстоянии е (эксцентриситет) относительно центра окружности С.. С.. 


Гипфциклойда Г о 
статора ротор 


Профиль 
статора 


- 150 -120 -90 60 90 120 150 
Рис. 1. Схема образования гипоциклоидального профиля торцевого сечения рабочего органа героторной машины, 
2, =4, К, =150 мм, К, =30 мм 
Внутри окружности С, перекатывается без скольжения окружность С; (производящая окружность 
гипоциклоиды ротора), одна точка которой образует гипоциклоиду ротора. Внутри окружности С, 
перекатывается без скольжения окружность С, (производящая окружность гипоциклоиды статора), одна точка 
которой образует гипоциклоиду статора. Заметим, что радиусы окружностей С; и С, равны эксцентриситету е. 


Профили ротора и статора образуются как эквидистантные линии, отложенные на расстояние К, от 
соответствующих гипоциклоид. 

В такой постановке профили ротора и статора обеспечивают идеальный контакт при эксцентричном 
вращении ротора. При необходимости обеспечить натяг/зазор в соединении следует изменить эквидистантное 


расстояние А, для профиля ротора на величину А. 
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Перейдем к математическому описанию методики образования профилей торцевого сечения ротора и 
статора. Число зубьев статора равно 


2) =2+1[. (1) 


Радиус основной окружности гипоциклоиды статора К, 


р5 


определим из условия (рис. 1): 
Ё,, =А, В (2) 


Определим эксцентриситет е передачи по условию равенства длины окружности С, целому числу 2, длин 
производящей окружности С, радиусом е: 
2рК 


 — 2 рег, › 
тогда 
е=К 


5 


И (3) 


Заметим, что при 2 =1 (механизм Муано) из условия размещения ротора в полости статора необходимо 


определять эксцентриситет и радиус ролика К, по другим формулам: 
е= К, /ЗиК, =2е, если 2, =1. (4) 


Радиус основной окружности гипоциклоиды ротора А,. равен (рис. 1) 


В. =А,-е. (5) 
Введем функцию для определения радиуса основной окружности ротора или статора, выполняющую 


условия (2) и (5): 


Бе 


К(2)=В, —В,-е.(2,-—2). (6) 
Параметрические формулы гипоциклоид ротора и статора в декартовой системе координат ХОУ: 


х(ф,2) = (=). | (1-1/2)соз(ф/=)+ (а/ 2) соз((1-1/=)$) | 


У(Ф,=) = В(2)-[ (1—1/ 2) (Ф/=)- (а/ 2) ((1-1/2)9) | | ы 


где 2 — число ветвей гипоциклоиды ротора 2, или статора 2,; а — коэффициент укорочения гипоциклоиды. В 
случае обыкновенной гипоциклоиды в теории а=1. Однако в практических расчетах с применением систем 
компьютерной математики (Ма сад, МАТНГАВ) при а=1 возникают сложности с вычислением 
эквидистантных профилей ротора и статора. Поэтому в модели мы допустили а = 0,99, что для достижения 
наших целей вполне оправданно; ф — угол поворота радиус-вектора, проведенного из начала координат до 
каждой точки гипоциклоиды. Если гипоциклоида должна быть построена по №, точкам, то 


Ф(1,2) =2-п:1:2/М,, 1=0...М,. (8) 


р? 


Профили ротора и статора определяются по известным формулам эквидистантной траектории относительно 
гипоциклоид ротора или статора на расстояние К, с учетом натяга в сопряжении А 


Хула [22 +92 
(9) 


оф) о- (В+). [2 +82 


При работе передачи центр ротора движется по окружности радиусом е и поворачивается на угол 9=0...2м. 
Все точки ротора при этом совершают дополнительный поворот на угол \у, определяемый функцией 


лу(9) = и. (1-2, /2)). (10) 


Модель вращения ротора описывается параметрическими функциями преобразования координат при 
перемещении и вращении объекта: 


Х, (0, ф, А) =е-с0$(0) + Х,(ф,2,А)соз(цу(9)) — 7,(ф,2,,А)зи(у(а)) 


У (а, Ф.А) =е-з(о)- (а, = 1,0, Х (6,2, А) зу (о), (2, А)соз(у(о))" вы 
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В выражении (11) применена условная функция, обнуляющая перемещение по осиУ в случае 
прямолинейного перемещения ротора в полости статора при <, =1. 

Опыт применения данной методики показал, что эквидистантные профили ротора и статора при увеличении 
эквидистантного расстояния или радиуса ролика К, могут потерять гладкость. Это происходит в точках 
перехода от эквидистанты окружности в эквидистанту гипоциклоиды. На рис. 2 представлена проблема потери 
гладкости линии профиля статора в окрестности +1 мм от точки перехода эквидистантных линий. 


В математике известна формула критерия гладкости функции, заданной в параметрическом виде. Для 
профиля статора такая функция будет иметь вид: 


_ 07, (р,г,0) ОХ „(р,т.0) 
в. бр 


Е(р.&) (12) 


) 


где & = <,; р — угол радиус-вектора точки профиля статора. 


Необходимо отметить, что функция РЁ, во всем диапазоне (1, 5.) ведет себя нестабильно и найти какие- 


5 


либо особенности ее поведения в интересующей нас точке перехода эквидистантных линий достаточно сложно. 
Поэтому поведение функции Е изучалось только в окрестности точки перехода при изменении углового 


5 


параметра 
р=п/ <, -т...п/ <, +т, ИЕ 
где т — угол отклонения от угла радиус-вектора точки перехода. 
т=2.-р-и, (14) 
где и — часть углового сектора профиля статора для изучения гладкости в окрестности точки перехода. 
Величину и (для <, = 2..12) определяем по эмпирической формуле: 
и=0,12 —0,01-(5, —1). (15) 


На рис. 2 можно заметить, что с увеличением радиуса ролика в окрестности точки перехода эквидистантных 
линий профиль статора из гладкого превращается в самопересекающийся. Поведение функции Ё также 


изменяется, и при появлении эффекта самопересечения профиля возникает петля (рис. 2 в). 

Этот эффект положен в основу предлагаемого авторами статьи алгоритма контроля гладкости профиля 
статора, приведенного на рис. 3. Следует отметить, что, судя по многочисленным замерам, если критерий 
гладкости статора выполняется, то гладкость профиля ротора проверять не требуется, так как он построен на 
гипоциклоиде с меньшим значением радиуса основной окружности. 


Ур, Ур, Ур, 
мм мм мм 
40.49 Точка перехода 50,15 Точка Е 59.80 Точка перехода 
39,15 48,81 58,46 
Ал 1117 112,3 Хр, мм 103,8 104,4 105,0 Хр, мм 96,67 97,34 98,01 Хр, мм 
Е, Е, Е, 
0,3 4 Гладкость есть 0,34 Гладкость еще есть 0,34 Гладкость нарушена 
| 033 А | 
= -1 -1 
1 1117 1123 Хр, мм 103,8 104,4 105,0 Хр, мм 96,67 97,34 98,01 Хр, мм 
а) 6) в) 


Рис. 2. Профиль статора в окрестности +1 мм от точки перехода эквидистантных линий (верхний график) и значения 
функции РЁ; (нижний график) при <, =З и К, =150 мм, для К, :а— 40мм; б — 50 мм; в — 60 мм 
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Ар(р,2,4),2,,т,М,,АФ 


Гладкость 
нарушена 


Гладкость не 
нарушена 


| Хр(р,2,,0) > 


Ар(р+АФ,2,,0) 


Рис. 3. Схема алгоритма контроля гладкости профиля статора 


Контроль гладкости профиля статора должен проводиться всегда при проектировании геометрии торцевого 
сечения героторной машины. С целью изучения параметра максимального радиуса ролика у ниже 


г тах > 
которого гладкость кривой не нарушается, авторы провели серию численных экспериментов, в которых 
задавались различными значениями контурного радиуса К, и числами зубьев ротора 2|. 

Первым целевым параметром данного исследования является площадь живого сечения рабочего органа 5. 
Этот параметр входит в формулу подачи героторной машины и определяет ее производительность. С 
имеющимися параметрическими функциями, описывающими профили ротора и статора, задача по определению 5 
сводится к вычислению разности площадей параметрических кривых статора 5, и ротора 5, (9). Формула 


площади параметрической кривой известна из математики, и поэтому в методике авторов она выглядит так: 


ОУ, (р,2,,0) ОУ, (р,2,А) р 
д 


222п 
$5=5,-5, =] Х,(р,2.,0)- бр 


ар- [| Х,(р,2,А). (16) 


Вторым целевым параметром исследования является приведенный радиус Ри кривизны при контактном 


взаимодействии ротора и статора. Напомним, что согласно теории Герца, при контакте двух цилиндрических 
поверхностей рассчитывается приведенный радиус, который В формуле контактных напряжений находится в 
знаменателе. Следовательно, величина контактных напряжений наряду с физико-механическими параметрами 
материалов контактируемых поверхностей будет иметь наибольшую величину при наименьшем значении 


приведенного радиуса р„,, определяемого по формуле: 


р„ = (р, -р,)/ (р, +р,), (17) 


где р, — радиус эквидистантной окружности ротора, равный параметру А,; р, — радиус кривой профиля 
статора. 
Заметим, если р, >|, то числитель в формуле (17) всегда растет быстрее, чем знаменатель, следовательно, 


наименьшее значение Рир ви Полностью определяется параметром наименьшего радиуса статора ри. 


В торцевом сечении рабочего органа героторной машины при любом угле поворота ротора осуществляется 
непрерывный контакт его эквидистантной окружности со сложной кривой профиля статора. По результатам 
наблюдений за динамикой контактных взаимодействий можно сделать вывод, что наибольшие напряжения 
возникают в точках перехода эквидистанты гипоциклоиды в эквидистанту окружности профиля статора. 
Следовательно, наименьший радиус линии профиля статора р. „также будет находиться в зоне перехода 


эквидистантных линий и определится по известной из математики формуле: 


3 
22 


аХ,(р,2,,0)\' [7,(р,2,,0) 


Ор бр 
Рыивы В) = , (18) 
ОХ (р,2,,0)\ [627,(р.2,,0)\_ [22Х „(р,2,,0)\ [0,(р,2,,0) 
Ор др? др? бр 
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где р=п/2, угловой параметр точки профиля статора в месте перехода эквидистантных линий. В случае 
потери гладкости линии статора параметр р., „„ стремится к нулю. 


Описанная выше методика реализована в системе математических расчетов Ма@са4. При выполнении 
расчетов коэффициент укорочения гипоциклоиды а=0,99, №, = 40 000. 


На первом этапе исследований планируется изучение максимального радиуса ролика у при котором 


г тах > 
критерий С’ гладкости линии профиля статора равен нулю. Изменяемые параметры — контурный радиус А, 
(100, 150, 200 мм) и число зубъев ротора2, от 1 до 7). Расчеты проводятся при последовательном увеличении 


фиксируем с 


параметра А, от 5 мм до срабатывания критерия гладкости С. Исследуемый параметр К, и. 
точностью до 1 мм. 


На втором этапе проводятся серии расчетов для изучения влияния радиуса ролика В, (от 5 мм до у 


числа зубьев ротора 2, (от 2 до 7) и контурного радиуса К, (100, 150, 200 мм) на площадь живого сечения 5 и 
минимальный приведенный радиус контакта рот. 


Описанная методика расчета геометрических параметров профиля торцевого сечения рабочего органа 
героторной машины может быть реализована в любой другой системе программирования или в системе 
компьютерной математики. 

Результаты исследования. Данные численных опытов по определению параметра максимального радиуса 
ролика ^ при котором гладкость линии профиля статора сохраняется, приведены в таблице 1. 


7 тах > 


Таблица 1 


Максимальный радиус ролика К (мм) при различных контурных радиусах К, (мм) 


’ тах 


2 К, =100 К, =150 В, =200 
1 66,6 100 133 
2 45 68 91 
3 27 41 55 
4 20 31 40 
5 15 23 31 
6 12 19 25 
7 10 16 21 


На основании данных таблицы 1 построены графики (рис. 4 а). Видно, что оцениваемый параметр К при 


г тах 
увеличении контурного радиуса № растет пропорционально. Введем в рассмотрение относительный 
параметр А, = А, „‚ /А,, график которого представлен на рис. 4 6. 


г тах 


г тах 2 мм К 
140 0,9 
120 0,8 
0,7 1 
100 р 
0,6 
80 3 0,5 
60 № 0,4 
2 
40 9 2 
1 0,2 
20 о, 
0 | 0 
1 3 5 7 2, шт 1 6 И 16 2 Ш 
а) 6) 


Рис. 4. Зависимость максимального радиуса ролика от числа зубьев ротора: 


а— В, „при К, : 1 — 100 мм; 2 — 150 мм; 3 — 200 мм; б — К, (1) илиния тренда (2) с экстраполяцией до 2, =16 


г тах 
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Эмпирическая формула, Сс помощью которой можно определить максимально допустимый радиус ролика 


у. имеет вид: 


7 тах > 


В, =0,7682,0°% (19) 


При 2= 15, А, = 150 мм, К, = 0,0517 мм, К = 7,76 мм гладкость профиля достигается при К,= 7,5 мм. 


г тах 


При 2 =10, К, =200 мм, КД, = 0,0775 мм, В 


т тах 


= 15,5 мм гладкость профиля достигается при К,= 15,5 мм. 


Таким образом, можно сделать вывод об удовлетворительной точности эмпирической формулы (19) для 
предсказания максимально допустимого параметра К, пах. 

Следует отметить, что формула (19) получена при значениях коэффициента укорочения а = 0,99. 
Уменьшение а приводит к достижению гладкости профиля и при более значительном увеличении 


параметра К,. Однако при этом изменяется форма зуба циклоидального зацепления и нарушаются условия 
контактного взаимодействия ротора и статора, особенно в тех местах, где ротор должен касаться статора 


боковыми сторонами зубьев. Расчеты показывают, что при а=0,99 и К, = 38 мм критерий гладкости 


находится на пределе выполнимости. Однако при а=0,99 и К, = 38 мм достигается хорошая гладкость 
профиля статора, но нарушается кинематический контакт деталей рабочего органа машины, восстановить 
который удается лишь специальным подбором значений параметров гипоциклоид ротора и статора. Поэтому 


исследование влияния коэффициента укорочения на параметр у выведены за рамки данного исследования. 


г тах 

В таблице 2 представлены Данные расчетов целевых параметров исследований для числа зубьев 
ротора 2 = 3 при различных значениях контурных радиусов. Полученные значения площади живого сечения 5 
и минимального приведенного радиуса контакта Рир п В Последних двух колонках таблицы 2 приведены к 


безразмерному виду для удобства анализа. 


Таблица 2 
Данные расчетов целевых параметров при 2!=3 
В, мм Риршы» ММ 5, мм? 5Лпах (5) Рурк ВХ (Руна) 
К, = 150 мм 
5 4,72 18038 0,9712 0,2656 
10 8,83 18256 0,983 0,4969 
15 12,28 18415 0,9915 0,6911 
20 14,98 18513 0,9968 0,843 
25 16,85 18572 1 0,9482 
30 17,77 18530 0,9977 1 
35 17,63 18450 0,9934 0,9921 
40 16,28 18307 0,9857 0,9162 
К, =200 мм 
10 9,14 32275 0,9786 0,3857 
20 16,37 32738 0,9926 0,6909 
30 21,36 32982 1 0,9015 
35 229 32980 0,9999 0,9664 
40 23,69 32943 0,9988 0,9998 
45 23,695 32844 0,9958 1 
50 22,82 32685 0,991 0,9631 
55 20,98 32466 0,9844 0,8854 
К, = 100 мм 
5 4,57 8068 0,9773 0,3857 
10 8,18 8184 0,9914 0,6903 
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В, ‚ ММ Рир ши 2 мм 5, мм? 5Итах (5) Рир ши ? пах ( Рир тт ) 
К, = 150 мм 
13 9,83 8225 0,9964 0,8295 
15 10,68 8245 0,9988 0,9013 
И 11,32 8255 1 0,9553 
20 11,85 8235 0,9976 1 
24 11,64 8189 0,992 0,9823 


На рис. 5 представлены графики изменения целевых параметров исследований при 2=3. На них нетрудно 
заметить, что при увеличении радиуса ролика А, происходит незначительное изменение площади живого 


сечения 5 и достаточно сильное изменение приведенного радиуса контакта, максимальное значение которого 


приходится на достаточно узкий диапазон К. , что представляет для проектировщика наибольший интерес, так 
как при этом контактные напряжения будут минимальны. Расчеты для других значений 2, (от 2 до 7) показали 


примерно такую же качественную картину, как и на рис. 5. 


1,0 "ТТ 1,0 »-- 97932 9:$° 1,0 о. 9 +: 90-0 
Ро р 0,9 И. 5 % | 05 % р 
0,8 |. 1 е 1 Г 
Ед 0,8 ; о ; 
0,7 р. а : 
р : 0,7 — 0,7 я 
и $ 0,6 } 0,6 
РИ ь 0,5 ; 0,5 
р 0,4 | 0,4 : 
0,3 | 5 ( 
0,2 0,3 0,3 
0 10 20 30 Юмм 0 20 40 Кьмм 0 10 20 КЮ мм 
а) 6) 6) 


Рис. 5. Зависимости целевых параметров исследований от радиуса ролика К, при 21 =З, а = 0,99, Д=0: 


1 — 5/пах(5); 2 — Рив /мах(Рир ви ) (таблица 2); К, ‚мм: а— 150; 6 — 200; в — 100 мм 


В таблице 3 представлены Данные расчетов оптимальных радиусов роликов Вор для различных 
параметров 2 и К, . В результате анализа выяснилось, что параметры Вор и К, для одного 2 изменяются 


практически одинаково. Поэтому, если ввести параметр относительного оптимального радиуса ролика 


в, =, /В,, (20) 


г ор! 


то с помощью однофакторной зависимости вида К, =/(2,) можно на предварительном этапе проектирования 
рабочего органа героторной машины определять оптимальный радиус ролика К... 

В результате аналогичных рассуждений для приведенных в таблице 3 параметров 5 и р,„„„„ введены также 
безразмерные параметры относительного приведенного радиуса контакта р, и относительной площади 
живого сечения 5, : 


ПЯТ, ре ы Г (21) 


приш 


Новые относительные величины также представлены в таблице 3: 


Киреев С.О. и др. Оптимизация геометрических характеристик циклоидальных профилей героторных гидромашин 


Таблица 3 
Оптимальные радиусы К,„„ и относительные параметры зацепления 
В,, мм 2, К, р» Мм | 5, мм? Ррши» ММ К, Рура 5, 
200 2 68 36103 43,2 0,340 0,216 0,903 
150 2 51 20318 32,2 0,340 0,215 0,903 
100 2 34 9015 21,65 0,340 0,217 0,902 
100 3 21 8235 11,85 0,210 0,119 0,824 
150 3 32 18530 17,8 0,213 0,119 0,824 
200 3 42 32800 23,1 0,210 0,116 0,820 
200 4 28 30015 15,32 0,140 0,077 0,750 
150 4 21 16880 11,5 0,140 0,077 0,750 
100 4 14 7495 7,63 0,140 0,076 0,750 
100 5 10 6860 5,35 0,100 0,054 0,686 
150 5 15 15441 8,1 0,100 0,054 0,686 
200 5 20 27452 10,81 0,100 0,054 0,686 
150 7 9 13050 4,7 0,060 0,031 0,580 
200 7 12 22100 6,2 0,060 0,031 0,582 
100 7 6 5780 3,2 0,060 0,032 0,578 


Следует отметить, что все вновь введенные относительные величины с высокой степенью достоверности 


подчиняются аппроксимирующим формулам, включающим один фактор — число зубьев ротора 2. 


Формула относительного оптимального радиуса ролика К, имеет вид: 


К,, =0,9331-2. "3 (2 =0,9965), 


где К? —_ величина достоверности аппроксимации. 


Формула относительного приведенного радиуса контакта р„,, имеет вид: 


Формула относительной площади живого сечения 5 имеет вид: 


бы = 0635 (В 0.9995) 


5, =10175—0,06412, (Е? =0,9911. 


(22) 


(23) 


(24) 


Рассмотрим пример применения полученных результатов при проектировании геометрии профиля 
торцевого сечения героторной машины. Зададимся следующим набором исходных данных: контурный 


радиус К, = 120 мм, число зубьев ротора 2, = 4, натяг в сопряжении Д=0 мм (эмпирические формулы получены 


при нулевом натяге). Относительный оптимальный радиус ролика А,, по формуле (22): 


Абсолютное значение радиуса ролика К, = А „„-К, = 0,135.120 =16,2мм. 


К, =0,9331-4-13° = 0,135. 


Относительный приведенный радиус контакта р„,, по формуле (23): 


р», = 0,635-4-*5%* = 0,076. 


Абсолютное значение минимального приведенного радиуса контакта: 


Относительная площадь живого сечения 5, по формуле (24): 


р пришш 


=’ 


про 


К, =0,076.120 =9,12 мм. 


5, =10175-0,0641.4 = 0,761. 
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Абсолютное значение площади живого сечения: 
$5=5,.В?= 0,761-120? = 10958 мм. 


Расчеты по программе Маса4 дают следующие величины: 5 =10820мм*, Р‚ршы = 9,167 мм. Как видно на 


этом примере, эмпирические формулы достаточно точно предсказывают основные параметры геометрии 
торцевого сечения героторной машины. 

Обсуждение и заключение. Проведенное исследование направлено на решение задачи совершенствования 
рабочих органов героторных гидромашин и в достаточной степени соответствует работам отечественных и 
зарубежных авторов [10-15]. Основное отличие заключается в том, что при построении профилей ротора и 
статора использованы идеальные гипоциклоиды, обеспечивающие наибольшую кинематическую точность 
движения и контактного взаимодействия элементов зацепления. Известную проблему обеспечения натяга в 
сопряжении предлагается решить за счет разности параметров радиуса ролика, а для контроля гладкости 
профилей разработан специальный алгоритм [12]. Оптимальную геометрию профиля, как и в работе [10], 
предлагается определять на основе максимальной площади живого сечения рабочего органа героторной 
машины. Но при этом приведенный радиус контакта также должен иметь наибольшую величину, как залог 
наименьших контактных напряжений. В данном исследовании на основе большого числа расчетов по 


специально разработанной программе Маса4 в диапазоне чисел зубьев ротора 2 от 2 до 7 и контурного 
диаметра К, от 100 до 200 мм были определены оптимальные значения эквидистантного расстояния 
гипоциклоид, или радиусов роликов К,, при которых достигается наилучшее сочетание приведенного радиуса 


контакта и площади живого сечения рабочего органа героторной машины. Выяснилось, что введение 
относительных величин радиуса ролика, площади живого сечения и приведенного радиуса контакта позволяет с 
помощью простых эмпирических зависимостей на раннем этапе проектирования получить абсолютные 
значения этих же величин для любого контурного радиуса в пределах установленных диапазонов. 

Разработанная методика дает возможность на этапе проектирования профилей ротора и статора оценить 
производительность и работоспособность героторной машины. В перспективе подобные изыскания могут быть 
проведены для эпициклоидных и смешанных зацеплений в героторных гидромашинах. 

Результаты исследований могут быть использованы в машиностроении при проектировании героторных 
гидравлических машин для улучшения их технико-эксплуатационных характеристик. 
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